Magnetic anisotropy of vicinal (001) fee Co films: role of crystal splitting and 
structure relaxation in step-decoration effect 
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The uniaxial in-plane magnetic anisotropy (UIP-MA) constant is calculated for a single step on 
the (001) surface of fee Co(A'') films. The calculations are done for both an undecorated step and the 
step decorated with one or more, up to 7, Cu wires. Our objective is to explain the mechanisms by 
which the decoration decreases the UIP-MA constant, which is the effect observed experimentally 
for ultrathin Co films deposited on vicinal (001) Cu surfaces and can lead to reorientation of mag- 
netization within the film plane. Theoretical calculations performed with a realistic tight-binding 
model show that the step decoration changes the UIP-MA constant significantly only if the splitting 
between the on-site energies of various d-orbitals is included for atoms located near the step edge. 
The local relaxation of atomic structure around the step is also shown to have a significant effect 
on the shift of the UIP-MA constant. The infiuence of these two relevant factors is analyzed further 
by examining individual contributions to the UIP-MA constant from atoms around the step. The 
magnitude of the obtained UIP-MA shift agrees well with experimental data. It is also found that 
an additional shift due to possible charge transfer between Cu and Co atoms is very small. 

PACS numbers: 75.30.Gw, 75.70.Cn, 75.70.Ak 
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I. INTRODUCTION 

The magnetic anisotropy (MA), determined by the de- 
pendence of energy on the orientation of magnetization, 
is one of the basic properties of magnetic systems. A 
non-zero MA energy arises due to two different physi- 
cal effects: long-range magnetic dipole-dipole interaction 
and spin-orbit (SO) coupling. The latter interaction cou- 
ples the spin of an electron to the local electric field as 
the electron travels through the system which makes the 
interaction sensitive to the local atomic structure. This 
leads to the magnetocrystalline (MCA) component of the 
MA energy. As the symmetry is reduced in nanoscopic 
structures, the MCA energy (per atom) becomes larger, 
normally by two orders of magnitude, than in bulk cu- 
bic crystals. ^■■^ An example of such low symmetry struc- 
tures is a ferromagnetic ultrathin film deposited on a vici- 
nal surface which consists of monoatomic steps separated 
by flat terracesiS. Vicinal surfaces can be obtained by 
slightly miscutting low-index [e.g., (001)] surfaces of sub- 
strate crystal. 

In the present paper, we study theoretically MA of fee 
Co films placed on a vicinal (001) Cu surface with step 
edges running along the direction [110]. In such systems, 
the magnetization M lies within the Co film plane due 
to the negative sign of the out-of-plane MA constant K2 
which determines the energy dependence K2 sin^ 9 on the 
magnetization polar angle 9; cf. Ref. 0- With M lying 
in-plane (i.e., for 9 — 7r/2), the system energy depends 
on the azimuthal angle </) which M makes with the [110] 
axis: 

E{(t>) ^Ki/As\r?{2(t>)+K^s\T?4>-MHcos{(t>-(t)H). (1) 



Here, the first term on the right-hand side is due to bi- 
axial (fourfold) MA which arises because the directions 
[110] and [100] are non-equivalent in systems with cubic 
crystal structure. The cubic MA constant Ki is positive 
for fee Co/Cu(001) filmsi^Thus, in flat Co/Cu(001) films 
the magnetization lies along one of the two equivalent 
directions: [110] or [110] , the easy axes, while the direc- 
tions [100], [010] are hard axes. The existence of steps in 
vicinal Co/Cu(001) films introduces an additional, uniax- 
ial (twofold), magnetic anisotropy within the film plane. 
The directions [110] and [110] seize to be equivalent: one 
of them remains the easy axis, the other becomes the in- 
termediate axisi^ The sign of the uniaxial in-plane MA 
(UIP-MA) constant Ky^, equal to 



= E{(^ = V2) - E{c^ - 0), 



(2) 



a,t H — Q, determines the easy direction of magnetization: 
M is parallel (perpendicular) to the step edges when K^^ 
is positive (negative). 

The expression includes, optionally, an external 
magnetic field H which is assumed to lie within the film 
plane at the azimuthal angle (j>H to the [110] axis. Ex- 
perimentally, by applying the field H perpendicularly to 
the easy axis, one can measure the so-called shift field 
and determine the UIP-MA constant as 



(3) 



here Ms is the saturation magnetizationi^ This relation 
holds as long as is much smaller than Ki which is 
true for vicinal Co/Cu(001) films.^ 

The magnetic anisotropy of thin ferromagnetic films 
deposited on non-magnetic vicinal surfaces has been 
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the subject of many experimental papers. Apart from 
Co/Cu(001) systems,^iSiL2iiaJii2*H the reported work 
also concerned Fe films on vicinal Ag, Au, W, Pd, Cr 
(001) substratesi^iiSiiSiiLiSJ^ and vicinal Ni/Cu(001) 
systems.^^ Related experimental measurements of UIP- 
MA have been performed for Co/Cu(110) ultrathin 
filmsSSiSiS4 and also for Co films obtained by deposit- 
ing Co atoms at an oblique angle on a flat (001) Cu sur- 
face which leads to some structural anisotropyiS^ In the 
experiments on MA in vicinal films, it was investigated 
how the shift field Hs depends on the ferromagnetic film 
thickness, the miscut angle a and the miscut directioni^ 
These dependencies were found to be determined by the 
kind of the deposited film, its lattice type, the choice of 
substrate, and, in consequence, also by the strain present 
in the film. For Co/Cu(001) systems, the UIP-MA con- 
stant has been found to depend linearly on the mis- 
cut angle a which implies that is proportional to the 
step density; cf. Eq. Q below. The experimental de- 
pendence on the Co film thickness d is more complicated: 
the usual breakdowni^ = K^ + K^/d into the volume 
and surface terms is not generally valid, presumably, due 
to the anisotropic relaxation of the biaxial strain present 
in Co/Cu systems. 

Covering a vicinal ferromagnetic film with non- 
magnetic adsorbates, like Cu, Ag, O or CO, usually has 
a large effect on UIP-MA. ^•^■^^ In particular, it has been 
found experimentally that even submonolayer amount 
of Cu deposited on a vicinal Co film can decrease 
so significantly that, in some cases, becomes neg- 
ative and the magnetization direction switches to per- 
pendicular to the steps i^iSiSiiSiiiiiS. This effect is usu- 
ally explained by change in local MCA due to Co-Cu 
orbital hybridization at Co step edges after they are dec- 
orated with Cu atomS:^^^ Such qualitative explanation 
seems plausible, however it still requires further inves- 
tigation since few, and only partly relevant, theoretical 
calculations have been done so far. In Ref. 0, the en- 
ergy of a single step on a (001) fee Co film with non- 
collinear magnetic structure was calculated within an ab- 
initio model without SO coupling included. However, 
no values for the UIP-MA constants nor definite conclu- 
sions about the magnetization reorientation were given 
there. The relativistic full-potential linear-augmented- 
plane-wave calculationspS performed for a planar net- 
work of parallel Co wires separated with empty wires 
(i.e., a striped monolayer), have shown that magnetiza- 
tion changes its orientation upon filling the empty wires 
with Cu. However, this structure seems too simplified 
to extend the obtained results to the case of Co vici- 
nal films. In Ref. 30, more realistic systems, i.e., un- 
supported stepped Co and Fe monolayers were consid- 
ered using a tight-binding (TB) model, but the effect of 
step decoration was not studied. The present authors 
reported^ a TB calculation for a monoatomic step on a 
5-monolayer (ML) (001) fee Co film. It was found that 
the MCA contribution K^'^^ to the UIP-MA constant 
i^u was only slightly shifted after step decoration and 



the shift was much smaller than experimental values (cf. 
Sec. Unil. With the same TB model, shifts of by 
similar or smaller magnitude were found for steps on Co 
films with thicknesses other than 5 ML (in the range 2 
up to 10 ML). 

The present theoretical study focuses on explaining 
the mechanisms by which decoration of Co steps with 
Cu atoms changes UIP-MA in vicinal Co/Cu(001) sys- 
tems. We calculate K^'^^ for stepped Co films using an 
improved TB model. This model includes crystal field 
splittings (CFS) introduced into the on-site energies of 
d-orbitals on every atom in the system. Similar, though 
simplified, CFS have previously been included and shown 
to play a crucial role in calculations of the out-of-plane 
MA constant of flat films and multilayers iSSiSi^iSi^SiSi 
Here, we also investigate how the UIP-MA constant 
is modified by the relaxation in the Co film structure 
around the step. In addition, a potentially significant 
effect of the charge transfer (CT) between Cu and Co 
atoms in the decorated steps is studied. We do not 
include the contribution to UIP-MA coming from the 
dipole-dipole interaction because it remains almost un- 
changed when the Co steps are decorated with non- 
magnetic atoms (cf. Sec. IIIII Ref. 38). 



II. MODEL 
A. Geometry. Structure relaxation 

In experiment, cobalt films deposited with molecu- 
lar beam epitaxy grow pseudomorphically on the vici- 
nal (001) Cu surface for film thicknesses larger than 2 
monolayers (ML) and less than 10-15 MLsi^. Thus, the 
Co film surface replicates the underlying stepped sub- 
strate surface quite well in this thickness regime. When 
submonolayer amount of copper is next deposited on 
such vicinal films, Cu atoms decorate the Co step edges 
as has been shown by scanning-tunneling-microscope 
measurements ^ 

The miscut angle a determines the terrace width and 
consequently the number iVrows of atomic rows (lying par- 
allel to [110]) on each terrace»i^ 

7V,ows = -1=^ (4) 

V 2 tan a 

The steps on vicinal surfaces are normally well sepa- 
rated. Indeed, Arrows is equal to approximately 400, 
25, 12 atomic rows for the vicinal Co/Cu(001) sys- 
tems with small miscut angles a = 0.1°, 1.6°, 3.4°, re- 
spectively, which were studied experimentally in Refs. 
5 6 1,9,10 11,12, Therefore, as the decoration of a given 
step is expected to change MCA only locally, we restrict 
our theoretical model to a single monoatomic step on Co 
(001) fee film. For further simplicity, we assume that 
the step is present only on the upper film surface and 
the bottom surface is flat so the system, hereafter called 
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FIG. 1: Geometric structure of (001) fee Co{N/N - 1) step: 
(a) undecorated, (b) decorated with one Cu atomic wire 
(shaded circle). 



Co{N/N — 1) step, is actually built of two joined semi- 
infinite slabs with thicknesses of N and A^— 1 monolayers, 
correspondingly [Fig. da)]. For the decorated systems, 
there are one or more Cu wires attached to the step edge 
on the upper surface on the Co film [Fig. Wi^)] ■ For later 
reference, we assume that the axes x, y, z are along the 
directions [110], [110] and [001], correspondingly, and the 
step-edge Co wire is located at the position a; = 0, z = 0. 

Due 1.8% mismatch between the lattice constants of 
cobalt and copper, the cobalt film is subject to biax- 
ial strain. This leads to tetragonal distortion of the Co 
film lattice: the in-plane distances between Co atoms 
are expanded, to match the Cu lattice constant, while 
the distances between Co planes become compressed in 
comparison to bulk fee Co lattice. In the present work, 
we consider Co{N/N — 1) steps of thicknesses < 10 
which have either (i) perfect or (ii) tetragonally distorted 
fee lattice. In the latter case, we allow for additional re- 
laxation of the atomic positions in the vicinity of the 
step. The relaxed atomic positions Ri are obtained by 
minimizing the system energy 



E{Ri,R2,...)^Eb{Ri,R2,---) + \ J2 Vp{R,-R,). 

(5) 

It includes the band energy and the pair potentials 
Vp which are given by the analytical formulae with pa- 
rameters fitted^ to reproduce various elastic properties 
either known from experiment for Cu or found with first- 
principle methods for Co and Co/Cu systems. This 
semi-empirical model, based on the second moment tight- 
binding approximation, was previously used to find strain 
relief and shape evolution for a Co island on Cu surface. 
Here, we use the model to calculate the displacements of 
atomic positions around the step edge. These displace- 
ments (magnified by a factor of 30) are depicted for the 
Co(4/3) step in Fig. |21 very similar shifts are found for 
Co(8/7) step. They are shown relative to the biaxially 
strained slabs, Co(iV) and Co(iV — 1), for which the in- 
terlayer spacings are also obtained with Eq. |SJ|. 
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FIG. 2: Relative displacements (arrows) of the atomic posi- 
tions in the biaxially strained lattice around the edge of the 
(001) fee Co(4/3) step: (a) undecorated, (b) decorated with 
one Cu wire. The arrows are 30 times longer than the actual 
displacements relative to the lattice. 

B. Tight-binding model. Crystal field splittings 

The electronic structure is described with an extended 
TB hamiltonian 



E T^^^.^^cl^^Cjau 



(6) 



Here, /i, v denote any of nine orbitals, of s, p, or d sym- 
metry, located on every atom (Co or Cu), labeled with 
index i or j; a is spin and c\ c are creation and annihi- 
lation operators, respectively. The hamiltonian depends 
on the magnetization direction {9, 0) through the SO in- 
teraction 



so[ 



(7) 



^i' aa' 



where the elements (/icr|i • S\va'^ axB the known analyt- 
ical functions of 0, <\) (cf. Ref. |4l|). We assume the SO 
constant = 0.085 eV for Co atoms^iS*^, and = 0.10 
eV for Cu atomsj^S-^ 

The elements T^^ describing electron hopping are 
calculated with the Slater-Koster (SK) formulae;^^ the 
sum over j in Eq. © includes the first and second near- 
est neighbors for every atom i. We assume the values of 
two-center hopping integrals found for paramagnetic fee 
bulk Co in Ref. 46 by fitting to the ah initio band struc- 
ture. When strain and structure relaxation are present, 
the hopping elements T"^ are obtained by using the 

canonical scaling law: \Ri — ilj ["''"'"' here I and V 
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are the orbital quantum numbers for orbitals fi and v 
centered on atoms i and j which are focated at the posi- 
tions Ri and, Rj , accordingly; cf. Ref. 113 The hopping 
two-center integrals between Co and Cu atoms are ap- 
proximated by the geometric means of the corresponding 
Co and Cu integrals. 

The generalized on-site energies V^^^ (also known as 
on-site potentials) 



v; 



IfJ.V 



CFS 



(8) 



consist of two terms. The first one, diagonal in fi, v, has 
a common value Uu-^a for all orbitals /i of given angular 
symmetry, s, p or d (corresponding to Z = 0, 1 or 2). The 
other term describes CFS, i.e., possibly different values 
of Vf^^ for various orbitals with the same /; in addition, 
non-zero elements V^^^ with n ^ may exist. The CFS 
arise because various orbitals centered on a given atom i 
are oriented differently with respect to neighboring atoms 
which contribute to the electron potential V°'{r) on and 
around the atom i. The TB parameterizations for bulk 
cubic crystals, like that in Ref. EH include CFS, as the 
splitting between the on-site energies of t2g and Cg or- 
bitals. However, CFS are often neglected in TB mod- 
els for s ystems with lower symmetry like surfaces (cf., 
e.g., Ref. l3Ci48|) even though the corresponding CFS are 
larger than in bulki^ The splitting between the on-site 
energies of d orbitals oriented out-of-plane and those ly- 
ing in-plane, located on atoms at the film surface (or at 
an interface) has previously been used in the MA cal- 
culations for thin films and multilayersi^^i'^^i^'^i'^^i^^i'^'' It 
has been shown that such approximate CFS have a large 
effect on the MA energies^SiS 

In the present paper, we calculate V^^^^ for d orbitals 
/z, V centered on an atom i as the sum^ 



V 



CFS 



j 



(9) 



over its first nearest neighbors; we neglect CFS for s and p 
orbitals and all sp, sd, pd off-diagonal terms V^^^^ . In Eq. 
©, the terms I'i^^y are given by expressions formally 
identical^ to the SK formulae where the two-center hop- 
ping integrals Vdda-, Vdd-K, VddS are replaced with similar 
parameters ^^^^^('ti, m'), ("i, m'), v'^J^{m,m') de- 
scribing CFS. These two-center CFS parameters depend 
on the types m and m! of the atoms i and j, correspond- 
ingly, each of which is either Co or Cu. 

This method of calculating CFS in TB models has 
been presented in Ref. |^ and used in some theoretical 
papersi^2iS^iS^iS A more general approach to the on-site 
energies is described in Ref. [H^ while a different method 
using the Fourier transform of the potential is applied 
to find the CFS in the bulk cubic crystals in an early 
work on the TB model by Callaway and Edwards The 
expressions for V^'^ and vf^^^j can be derived similarly 
as the SK formulae, i.e., by approximating the potential 
V'^{y) with the sum of spherically symmetric atomic po- 
tentials wj(|r — Rjl) and expressing the orbitals ji and 



u centered on atom i in terms of the orbitals with the 
angular momentum quantized along the axis joining the 
atoms i and j. Then, the on-site energies Vj^^, found 
as the matrix elements of the hamiltonian ^V^ -I- F°'(r), 
include the terms (j/i|wj|w) which, after neglecting their 
spin-dependence, are denoted as vf^^j in Eq. 0. One 
concludes here that the elements vf^^^^ are, in general, 
not symmetrical in i, j, and nor do the two-center CFS 
parameters, like v'^J^{m, m'), need to have the same val- 
ues for {m,m') = (Co, Cu) and {m,m') = (Cu, Co). 
Because, the term V^^^^ becomes V^^^^ + cd^i, under 



the uniform shift v^J^ 



CPS 



c {j] = a, TT, S), we can 



include v^Jf into the diagonal on-site term Uu-a, Eq. 
|S1 and assume that v^J^ = hereafter. We determine 
the values of v^J^ and v^J^ by fitting the energies ob- 
tained with our tight-binding model for paramagnetic Co 
and Cu monolayers to the ab initio bandsi^ As a result, 
we have found the values i'2j^^(™, ™) = —0.17 eV and 
v^J^{m,m) = —0.10 eV, valid for both m = Co and 
m = Cu within the accuracy of the fitting method. The 
same values are assumed for the parameters v^J^{m, m'), 
v'^d^i'm^m') used to find the contribution of a Cu atom 
to the on-site energy V^^^ on a neighboring Co atom, 
i.e., for (m,m') = (Co, Cu), or vice versa, i.e., for 
(to, to') = (Cu, Co). Thus, we can obtain, e.g., the split- 
ting of 0.10 eV between the on-site energies of zx and yz 
orbitals for atoms located on the step edge. Similarly, we 
find that the orbitals 3z^ — and x^ — y^ are split by 0.20 
eV at the flat film surface; this value is close to the split- 
ting of 0.22 eV between out-of-plane and in-plane orbitals 
assumed previously for Co slabsi^^s^i We do not scale the 
two-center CFS parameters when the inter-atomic dis- 
tances change slightly in systems with relaxed structure; 
the canonical scaling law used for the hopping integrals 
seems not to be valid in this case. 

The on-site energies of down- and up-spin d orbitals 
on Co atoms are split by the exchange interaction: 



(0) 



(0) 



CX ' 



[l = 2). The exchange splitting 

A{^) ^ Abulk 
ex ex 71 ,f 



Ik 



(10a) 
(10b) 

(11) 



on a Co atom i is proportional to its spin moment Mi. 
Here, we assume that the exchange splitting in bulk fee 



Co is A 



bulk 



1.8 eV (cf. Ref. 



which, in the 



present TB model, leads to the bulk spin magnetic mo- 
ment Mbuik — 1-57 ^B, close to the ab initio values iWi^^ 
No exchange splitting is assumed for s and p orbitals on 
Co atoms and for all orbitals on Cu atoms. 

To find the yet undetermined on-site energy terms U^P 
we require charge neutrality, assumed separately for the 
d orbitals, 



bulk:; 



(12) 
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and the s and p orbitalSfSS 



The quantities n 



^bulk:m 



and n^^^ are the projected s, p, 



^bulk:m 



(13) 



and d occupations on an atom i while 
and nj^""^''" are the corresponding values found with the 
present TB model for the respective bulk metal m: fer- 
romagnetic fee Co or paramagnetic Cu. The conditions 
p2|l . are well satisfied at transition-metal surfaces 
as found in the ab initio calculations*'^ and have been 
used previously in TB models (cf., e.g., Refs. I3ll37ll6^ . 

According to Eqs. ((T^ . l(T^ . we usually assume that 
there is no CT in the investigated systems. However, for 
selected decorated Co(iV/A^— 1) steps, we study the effect 
of CT between Cu and Co atoms which leads to non-zero 
charges — (7i|e| on atoms i around the step edge. For 
this purpose, it is assumed here that a Co atom acquires 
qti- electrons from each of its first-nearest Cu neighbors 
so that — (7tr|e| is the charge transferred per one Cu-Co 
bond. In consequence, the corresponding electron occu- 
pations n^j^^ and nip , originally given by the right-hand 
sides of Eqs. ifT^ . ((T^ . have to be modified slightly. 

(i) 

So, we assume that only the d-orbital occupation n^^ 
changes, by the amount qi > 0, on Co atoms, due to 
a relatively large value of the local d-orbital-projected 
density of states (DOS) at the Fermi level. Further, it 



is assumed that the changes of the occupations 



and 



nip on Cu atoms are equal: An^'"" = A 



q^/2. It 



should be noted here that the actual distribution be- 
tween the sp and d orbitals of the electron deficiency 
qi = — 5gtr on the decorating Cu atoms is found to be 
irrelevant as far as UIP-MA is concerned. The approx- 
imate value qti- — 0.025 is obtained by fitting, within 
the applied TB model, the decrease Mi — Mbuik of the 
Co magnetic moment Mi at the Co/Cu interface in 
the (001) fee Co(6)/Cu(10) slabs to the ah initio value 
Mi - Mbuik = -0.06/^B found in Ref. IH 



Calculation methods: on-site energies, 
magnetocrystalline anisotropy 



The orbital-projected occupations n 



as 

well as the atomic moments M^, can be obtained easily 
once the diagonal part of the system's Green function 
(GF) G(z) = (z - iJxe)"^ is known,^lSi We determine 
the GF elements {ip,a\G{z)\i\xa) for a number of atoms i 
around the step by using the recursion method (RM).^"' 
We do this for all atomic wires which are distant by less 
than 1.5 lattice constants (ao) in the x ( i.e., [110]) direc- 
tion from the step edge; in each wire, we choose one atom 
with the smallest y (=0 or ao/(2-\/2)) coordinate. The 
on-site energy terms are then found for the chosen 
atoms by solving Eqs. ifT^ . ((T^ (appropriately modified 
for qtr 7^ 0) iteratively. For every iteration, the new value 
of U':^ for an atom i (as well as for all other, equivalent, 



atoms belonging to the same wire as the atom ?) is ex- 
pressed only with the values of and n'^ found in 
previous iteration steps. This means that, for each atom 
I, Eq. H12|l is solved as if it were a non- linear equation for 

a single variable . Similar assumption is made for Eq. 



(fT!^ and its dependence on J/^q and C/iV .^^ Despite such 
approximate way of solving the multi-atom set of Eqs. 
(|12(l , H13() , which are obviously inter- independent to some 
extent, only around 15 iterations are usually needed to 
fulfill Eqs. (|12|l . H13|l for all atoms i with the absolute er- 
ror less than 0.0001. A similar technique was previously 
successful for determination of layer on-site energies in 
thin magnetic films and multilayers The success of 
such algorithm proves a posteriori that the occupation 
tl"^ on an atom i depends mainly on the term J/j^^"* (i.e., 
for I — 2) determining the on-site energies of rf-orbitals 



on the same atom i; similarly, n. 



depends mainly 
onU^°^ andU^°^ (1 = 0,1). 

The occupations ni*' , rip , n|^*^ and the moments Mi for 
atoms i, which lie farther than 1.5ao from the step edge, 
are to a good approximation not affected by its presence. 
Therefore, the on-site energy terms U^^^ for these atoms 
are the same as for flat slabs TV or iV — 1 monolayers 
thick, and they are determined separately, prior to the 
main calculation for the Co{N/N — 1) step. 

The magnetocrystalline contribution to the UIP-MA 
constant K^, Eq. ||2Jl, is calculated with use of the force 
theoremSSiSi: 



2)- 



The electronic grand thermodynamic potential for a 
given magnetization direction {0,(j)), 



n{9, (f>) = -kBT J dE In 



1 + cxp 



kBT 



n{E), (15) 



is expressed in terms of the GF via the total density of 
states (DOS) 



n{E) = --Im Tr G{E + iS) 



(16) 



((5 — *■ 0-I-); ep is the Fermi energy and /cb the Boltzmann 
constant (in Eqs. ((T7|l exclusively, i denotes V^). 

We use thermodynamic potential Q for a finite temper- 
ature T (=300 K) instead of zero-temperature energy E 
because such approach has previously been found to help 
MA energies converge for flat thin films due to the ther- 
mal smearing of the Fermi level. "^'^'^^ However, it should 
be stressed that finite T is used only as a convenient tool 
in the numerical calculations. In particular, the present 
model does not account for the decrease of the spon- 
taneous magnetization which occurs in real systems at 
finite temperatures due magnetic excitations (cf. Refs. 
I3l37tl . 

We calculate the GF and hence the uniaxial in-plane 
magnetocrystalline anisotropy (UIP-MCA) constant us- 
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ing two distinct methods. The first one apphes the stan- 
dard RM to find the diagonal GF elements entering the 
expression for the DOS: 



n{E) 



--y2lm{jK\G{E + t6)\jK) 



(17) 



which comes from the general formula (|16|l when the 
trace is calculated with the states 



1 



V2 



(Ij/i T)±bM i». 



(18) 



These states are chosen instead of the states X)-, 
(bY* i)j because the continued fractions that express the 
diagonal GF elements in the RM converge much quicker 
in this case (cf. Ref. l68j) . However, still as many as 400 
levels of the continued fraction for d orbitals and 120 for 
s and p ones are required to calculate K^'^^ with the 
accuracy of 0.01 meV/(step atom) (i.e., per one atom 
on the step edgeU^ similar number of levels have been 
applied in Ref. '7(?. 

The UIP-MCA constant determined with RM is given 
by the sum of the atomic contributions 



MCA 



(i) 

MCA- 



(19) 



To reach the accuracy 0.01 meV/(step atom) for 



MCA 



it is sufficient to include the contributions fc. 



0) 



from 



^u;MCA 

all atomic wires which are within the 3.5 ao distance from 
the step edge (cf. Figs. El El [3 and the discussion below); 
again we choose one atom with the smallest y coordinate 
in each atomic wire. 

A many-fold reduction of the computational cost in 
the RM calculations of K^'^^, as weh as ni\ n^p \ n^'^ 
is obtained by truncating the Co(iV — 1/A^) step, nomi- 
nally infinite in-plane, to a finite cluster with the square 
base which is centered at x = y — and has its two sides 
parallel to the step edge. It has been found that the ef- 
fect of the introduced cluster boundaries is so small that 
j^MCA (jQgg jj^Q^ change by more than 0.01 meV/(step 
atom) provided that the cluster size in the x and y direc- 
tions exceeds 50 ao- This theoretical finding is interest- 
ing because the coefficients in the continued fraction are 
obtained through iterative application of the TB hamil- 
tonian i?TB, Eq. lO, on an initial state, like located 
on an atom close to the cluster center. Thus, these co- 
efficients start to be disturbed by the presence of the 
boundaries already after around 25 iterations for the as- 
sumed cluster size. Similarly, a cluster length and width 
of more than 30 ao guarantees that the errors in the occu- 
pations ni'^ , rip ^ , n^''' and the moments Mi are less than 
0.0001 for atoms lying close to the cluster center. This 
being true despite the fact that the number of levels in 
the continued fraction required for such high accuracy is 
30 for s and p orbitals, and 120 for d ones. 

Because, a few hundred levels of the continued frac- 
tion are needed to obtain well-converged results for the 



UIP-MCA constant, this calculation becomes difficult for 
Co{N/N — 1) steps with thicknesses N larger than 6 ML 
due to the long computation times required. To remedy 
this situation, a new method, based on the Mobius trans- 
formation (MT) approach^^ to calculate the GF, has been 
developed. First, the energy E (or grand potential fl) is 
shown to have the following general form: 

E{e, <p) - i;co„st + Sbuiki + E,,,{L). (20) 

where L is the length of the Go{N/N— 1) step after trun- 
cating it only in the x direction; the truncated system re- 
mains infinite in the y direction. The MT method allows 
us to determine each energy term, present in Eq. H20|l . 
separately. The term Eosc{L), is an oscillatory function 
of L and comes from the interference between the bound- 
aries of the truncated step. This goes to for the infinite 
step, i.e., for L oo. The bulk term i?buik(^, 0) is the 
sum of the energies of the disconnected flat slabs, N- and 
{N — l)-monolayers thick, and as such it contributes only 
to the cubic MA constant Ki in the in-plane dependence 
of the step energy on the magnetization direction {0, (j)), 
Eq. The constant term -Econst includes the contribu- 
tions -Eonds coming the boundaries of the truncated step 
(which are located a.t x ^ —L/2 and x ~ L/2). These, 
unwanted contributions can be found separately (by con- 
sidering truncated flat slabs), and thus eliminated from 
-Econst 7 leaving the the contribution Estcp coming the step 
region alone: 



E. 



step 



Ec, 



E, 



ends- 



Finally, the UIP-MCA constant is calculated as 



MCA 



E. 



'step V 



-)-E,top(<^=0) 



(21) 



(22) 



{6 = ^). More details of this MT approach can be 
found in Ref. |^ while the full description will be pub- 
lished elsewherei^^ The MT method is much more effi- 
cient computationally than RM and allowed us to ob- 
tained accurate values of the UIP-MCA constants for 
steps Co{N/N - 1) with thicknesses N up to 10 ML. 



III. RESULTS. COMPARISON WITH 
EXPERIMENT 

The atomic spin magnetic moments obtained within 
the described TB model by using RM are presented in 
Fig. 131 The moments at the Co film surface are en- 
hanced, in comparison to bulk, and similar results have 
been found in ab initio calculations. '^^i^^ The moment 
of Co atoms at the very step edge is further enhanced. 
This confirms the general rule that the local magnetic 
moment, as well as the local exchange splitting, are the 
larger the smaller the atom's coordination number, which 
is a consequence of the reduced local d-orbital-projected 
density of states due to the decreased hopping. A closer 
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FIG. 3: Atomic magnetic moments Mi (in (ib) for the (001) 
fee Co(4/3) step: (a) undecorated, (b) decorated with one Cu 
wire. The results are obtained within the TB model including 
neit/ier CPS nor structure relaxation (cf. Sec. IllBl . Selected 
Co atoms are marked with Roman numbers for reference in 
text. 

comparison of Figs. Ota) andl^Ib) reveals that the dec- 
orating Cu atom modifies the moments of the neighbor- 
ing Co atoms, due to the Co-Cu hybridization, as if it 
were a Co atom; at the same time the Cu atom acquires 
only a very small, induced, magnetic moment. In par- 
ticular, the moments of Co atoms labeled II and III in 
Fig. O become, after decoration, very close to the mo- 
ments of the atoms V and II, respectively, before deco- 
ration. This resembles the situation at the interface of 
Co and Cu films where the Co moments are very close to 
the bulk Co moments i^^'^^ Similar behavior and values 
of the moments around the step edge have been previ- 
ously found with tight-binding linear-muffin-tin-orbital 
method in Ref. 23- Presently, we have also found that 
the inclusion of CFS, biaxial strain and local structure 
relaxation, discussed above, modifies magnetization very 
weakly: atomic moments change by no more than 0.02 /xb 
in comparison to those shown in Fig. |31 

Slightly larger changes of Co atomic moments are in- 
duced by the CT between Cu and Co atoms if we as- 
sume that it takes place. Indeed, we find for the deco- 
rated Co(4/3) step that the moments on Co atoms I, II, 
III decrease from 1.821 /ib, 1-599 ^b, 1-722 ^b, respec- 
tively, for (jtr — 0, to the corresponding values 1-794 /iB, 
1-544 /iB, 1-668 /iB found with the fitted value gtr — 0-025 
of CT per Cu-Co bond; all these moments are obtained 
within the TB model with CFS and structure relaxation 
included. The moment change AMi induced by CT is 
very close to the negative value of the local electron ex- 
cess Qi on Co atom i, this being qi — qn- for atom I, 



0.6 - 




2 4 6 8 10 
N (ML) 

FIG. 4: UIP-MCA constant K^'^^ [cf. Eq. idJ] vs Co film 
thickness TV calculated for (001) fee Co{N/N — 1) steps: un- 
decorated (solid line) and decorated with one Cu wire (bro- 
ken line); cf. Fig. 1. The results are obtained within the TB 
model (cf. Sec. Ill Bl which includes: (a) neither CFS nor 
structure relaxation, (b) CFS but no structure relaxation, (c) 
both CFS and structure relaxation. 



and qi — 2qtr for atoms II, III. This happens because the 
condition An^j*' — qi has been assumed to hold for Co 
atoms, on which the minority-spin component of DOS 
at the Fermi level is much larger than the majority-spin 
onc- 

Because, the decoration of Co{N/N — 1) step with Cu 
atoms modifies the Co magnetic moments only in the 
immediate vicinity of the Co step edge (Fig. |3J), the re- 
sulting change in the dipole-dipole energy is minutei^ 
Thus, we can neglect the dipolar contribution isT^'P to 
the UIP-MA constant Ku in the present study which is 
focussed on the effect of the step decoration. However, 
in general, this contribution is not negligible: it is posi- 
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TABLE I: UIP-MCA constant K^^^ [in meV/(step atom)] 
of (001) fee Co(4/3) and Co(6/5) steps decorated with one 
Cu wire. The results obtained within the TB model which 
includes CPS, structure relaxation, and CT of gtr electrons 
per Cu-Co bond. 





qti = 


gtr = 0.025 


gtr = 0.05 


Co(4/3) 


0.18 


0.17 


0.15 


Co(6/5) 


0.33 


0.32 


0.28 



live and favors a magnetization lying parallel to the step 
edge. The calculations of K^^p for magnetic films on vic- 
inal surfaces and networks of magnetic stripes deposited 
on flat or vicinal surfaces will be presented elsewhere 

The calculations of the UIP-MCA constant K^'^^ 
have been done in three consecutive stages. The CFS are 
included in the TB model in the second stage and struc- 
ture relaxations in the third, while the initial TB model, 
i.e., in the first stage, includes neither of the two above 
ingredients. However, at every stage, the on-site energy 
terms Uu are found as described in Sec. El All the results 
for JC^^^ are gathered in Fig. 21 In each of the three 
variants of the TB model, the obtained K^"^^ is positive 
for both decorated and undecorated Co{N/N — 1) steps 
with thicknesses > 4. When no CFS nor structure re- 
laxation are present the step decoration with Cu has only 
small effect on K^'^^: it is changed by no more than 
0.05 meV/(step atom) and the shift is positivei^i The 
presence of CFS shifts the UIP-MCA constant if^^A 
the undecorated steps upwards by around 0.2 meV/(step 
atom); the inclusion of structure relaxation results in a 
further, slightly smaller, shift. At the same time, the 
average value of K^^^ for decorated steps remains al- 
most unchanged, being close to 0.3 meV/(step atom) for 
N > 6 ML, regardless of whether the CFS or structure 
relaxation are included. As a net result, the decora- 
tion decreases K^"^^ by 0.15 - 0.20 meV/(step atom) 
or 0.25 — 0.30 meV/(step atom) when only CFS or both 
CFS and relaxation, correspondingly, are included. 

Another mechanism which, in principle, could affect 
the UIP-MCA constant is CT between the decorating 
Cu atoms and the neighboring Co atoms. To study this 
possibility, we have calculated the UIP-MCA constant for 
the decorated Co(4/3) and Co(6/5) steps using the fitted 
value qti- = 0.025 and its double, qtr = 0.05. The results 
shown in Table |l] prove that CT has very small effect on 
the UIP-MCA constant for gtr = 0.025. For qtr = 0.05, 
the change of K^'^-^ is somewhat larger but it is still 
much smaller than the shift of K^'^^ due to CFS or 
structure relaxation. Thus, we conclude that CT is not 
responsible for the UIP-MCA shift induced by the step 
decoration and it can be neglected in the further discus- 
sion. However, this conclusion should not be generalized 
to other systems where larger interatomic charge trans- 
fers can take place and lead to more significant changes 
of MCA energies as a result of larger shifts of the local 
DOS with respect to the Fermi level. 



The obtained K^'^-^, shown in Fig. 0] oscillates with 
a period close to 2 ML when the step thickness N in- 
creases. Similar 2-ML oscillations of out-of-plane MCA 
constant have previously been found for flat Co films^ 
and Co/Cu multilayers.^^ It has been shown"^^ that these 
oscillations originate from the quantum well (QW) states 
present in the Co layer which have fe-vectors in the cen- 
ter of the two-dimensional Brillouin zone and energies 
close to the Fermi level. Similar QW states are likely 
to be responsible for the oscillations of K^'-"^{N) in the 
Co{N/N — 1) steps. Interestingly, such oscillations are 
absent once the atomic structure is allowed to relax (cf. 
Fig. EJc)); explanation of this behavior needs further 
study, especially, in connection to the absence of such 
oscillations in experiment iflii^ 

To compare the obtained results with the experimental 
data found for Co films on vicinal (001) Cu surfaces, we 
convert the measured Hg field, given in kA/m, into the 
UIP-MA constant Ku- 

Ku = 1.254 • IQ-'^NNrowsHs [meV/(step atom)]; (23) 

here N^ows is the terrace width which depends on the 
miscut angle a through Eq. Q. In Ref. 0, the Hg 
field is found experimentally for a — 3.4° and the Co 
film thicknesses N = 9.4, 10.8, and 13.2 ML. By tak- 
ing the corresponding values, AiJg = —13.5, —11.5, and 
-12.5 kA/m, of the differences AH^ = H^'''' - H^"'^'"' 
between the minimum Hg field found for the decorated 
films, and the field H™'^'"= for the undecorated 

film, we obtain with Eq. (|23() the respective experimental 
shifts AKu = -0.19, -0.19, and -0.25 meV/(step atom) 
of the UIP-MA constant. A very similar value of AK^ 
is found from the data reported in Ref. for 8-ML Co 
film deposited on slightly curved Cu surface which has 
a variable miscut angle a (< 5.4°). It was shown there 
that Hs depends linearly on a and it vanishes after de- 
position of an amount of Cu sufficient to decorate 70% 
of Co steps with single Cu wires. This result allows us to 
write down the equation: 

0.3i7r'*'=' + O.IH^'"' ^ (24) 

from which we find that H^"'' = -0A3H^'"^'"= and AH^ = 
— 1A3H^'^'^°'^ . Thus, from the experimental valuei^ 
^undoc ^ oc ^ 11.94 kA/m found for a = 3°, 
we get Ai^u — —0.23 meV/(step atom). The experi- 
mental values of the UIP-MA shift AK^ derived above 
are close to our theoretical results obtained both with 
or without the structure relaxation, provided that CFS 
are included. We regard the remaining discrepancy of 
0.05-0.01 mcV/(step atom) as very low since calculated 
MA energies, even those found with ab initio methods 
rarely match the experimental data perfectly. Good the- 
oretical prediction for AK^ confirms the conclusion of 
the experimental work reported in Ref. Ts' that the de- 
crease of Ku is 'caused only by the Cu adsorbates located 
near the step edges'. However, we do not obtain the sign 
change of after decoration which is usually observed 
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FIG. 5: Atomic contributions k'^^^MCA (bars) to UIP-MCA 
constant K^'^^ [cf. Eq. CHI] for (001) fee Co(6/5) step: un- 
deeorated (white bars) and decorated (grey bars) with one Cu 
wire. The horizontal dotted lines mark (001) atomic planes; 
atomic wires (Co or Cu) are located at points where the bars 
touch the lines which also serve as the zero reference levels for 
the bars. The plotted data are obtained within the TB model 
which includes neither CPS nor structure relaxation (cf. Sec. 



experimentall}i2i2iiiii^ and leads to switching of the mag- 
netization direction within the plane. Indeed, after the 
decoration, K^"^^ remains positive (for > 4); adding 
the dipolar contribution K^^^ will push the resulting 
even higher. This implies that the theoretical energies 
j{MCA g^j.g large in comparison with experiment, both 
for undecorated and decorated steps. This discrepancy 
is presumably due to some unaccounted terms of volume 
type which are clearly present in the experimental Hs {N) 
dependenciesi^ The origin of such terms is not cleariiiii 
and needs further investigation. It should be noted here 
that the tetragonal distortion of the biaxially strained 
fee Co film lattice, which is included in our calculations, 
retains the films fourfold in-plane symmetry and, thus, 
gives no direct volume contribution to UIP-MA>i2iii 

To develop a better understanding the origin of the 
obtained shift Ai^'^^^, we examine the contributions 
^u*MCA UIP-MCA constant from individual atoms 

around the step; see Figs. |31 El E| The largest contri- 

(i) 

bution ^u-MCA comes from the Co atomic wire at the 
very step edge, marked as atom I in Fig. |21 Other 
atomic wires that lie within the distance of 1.5 ao 
from the step edge, give smaller, but still quite sizeable 
contributions k^-MCA- T^^^ contributions from farther ly- 
ing atomic wires decay very quickly with increasing dw 
However, as already mentioned above, the contributions 
from several tens of atomic wires with < < 3.5ao 
are required to obtain if^^^ with the accuracy of 0.01 
meV/(step atom). The step decoration with Cu atoms 
affects UIP-MCA in three different ways. First, the ad- 
ditional electron hopping between the Cu atoms and the 



FIG. 6: Atomic UIP-MCA contributions k^^ljiQ^ obtained for 
(001) fee Co(4/3) step within the TB model which includes 
CPS but no structure relaxation (cf. Sec. I11B|I . Other details 
as in Fig. |3 
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FIG. 7: Atomic UIP-MCA contributions kl^-^cA obtained for 
(001) fee Co(4/3) step within the TB model which includes 
both CPS and structure relaxation (cf. See. Ill B|l . Other 
details as in Fig. 



neighboring Co atoms affects contributions around the 
step edge, mainly on atoms II and IV, for which the re- 
spective contributions ^u^mca shifted downwards af- 
ter decoration. However, shifts A/c^*j^(-,^ of both signs 
are found for several atoms around the step and, thus, 
they tend to cancel out to a large extent which explains 
why the resulting total shift Ai^^^^^, shown in Fig. EJa), 
is relatively small. The inclusion of CFS leads to a large 
shift of fcu-MCA on the step-edge atom (atom I), for the 
undecorated steps, while the contributions from other 
atoms remain almost unchanged. This means that the 
UIP-MCA is affected only very locally by CFS. Indeed, 
the effect of CFS on atom I is canceled once the step is 
decorated with Cu atoms; cf. Fig. El This happens be- 



10 



E 
o 

-t— • 
CO 

Q. 

B 

CO 

CD 

E 



< 
o 




CFS 



FIG. 8: UIP-MCA constant K^'^^ for the (001) fee Co(4/3) 
(solid eircles) and Co(6/5) (open cireles) steps decorated with 
one Cu wire; cf. Fig. mb). The results are plotted versus the 
multiplier Acfs that modifies the CFS contributions from the 
Cu atoms; cf. Eq. 1251 : see text. The results are obtained 
within the TB model including both CFS and structure relax- 
ation (cf. Sec. IllBt . The vertical dashed line denotes the 
nominal value Acfs = 1. 



cause the decorating Cu atoms remove, via Eq. (O, the 
splitting between the energies of in-plane orbitals, like yz 
and zx, located on atom I. On the other hand, the struc- 
tural relaxation has a more spatially extended effect on 
UIP-MCA: it results in small changes of the contributions 

^u;MCA 



for several atoms, both for undecorated and dec- 
orated steps. This is presumably the indirect effect of the 
tetragonal lattice distortion which changes the electronic 
structure. Furthermore, the local change of the geomet- 
ric structure induced by the step decoration modifies the 
shifts Afc^'j^Q^ also only locally, mainly on atoms I to V 
around the step. In each variant of the TB model, the 
total UIP-MCA shift, Aif^CA, can be found [with 0.01 

me V/ (step atom) accuracy] by adding shifts of ^u^mca 
from atoms with < dyj < 2ao; cf. Ref. |^ It sllould 
also be noted that the contribution to K^'^^ from the 
decorating Cu atom is negligible. 

In order to better understand the mechanism by which 
CFS affect UIP-MCA, we modify the contribution of 
the decorating Cu atoms to CFS on the neighboring Co 
atoms. It is done by multiplying the CFS two-center in- 
tegrals for m =Co, to' =Cu: 



Vddr, ("^, m') Acfs v^^^ (to, to') 



(25) 



with the common factor Acfs for i] — ct, 7t,6; cf. Ref. 

In particular, for Acfs — 0, when the Cu atoms do 
not contribute to CFS on the Co atoms, we find that the 
constant K^'^^ is larger by around 0.2 meV/ (step atom) 
than for the nominal value Acfs = 1 when the decorated 
steps Co(4/3) and Co(6/5) with relaxed structure are 
considered; cf. Figs. ElEIc). In fact, the value of K^'^^ 
for Acfs = becomes closer to K^'^^ obtained for the 
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FIG. 9: UIP-MCA constant K^'^^ for the (001) fee Co(4/3) 
step decorated with wcu Cu wires; the data point for wcu = 
corresponds to the undecorated step. The results are obtained 
within the TB model including both CFS and structure relax- 
ation (cf. Sec, nra . 



undecorated step in the TB model including CFS ; the 
remaining difference [0.08 meV/(step atom) for Co(6/5) 
step] is due to the change in the local structure relax- 
ation after the step decoration. We also note that with 
increasing Acfs , the UIP-MCA constant of the decorated 
steps decreases, almost linearly, however, it remains pos- 
itive for Acfs ^ 2. Thus, we do not consider possible 
inaccuracy in the derived two-center CFS parameters for 
Co and Cu as the source of the lack of magnetization 
switching in the present theoretical calculations. 

Finally, we study how the UIP-MCA changes when 
more than one, wcm Cu wires are attached to the Co 
step edge, in a similar way to that depicted in Fig. ^ 
The results obtained for < wqu < 7 Cu wires. Fig. 
13 confirm the experimental finding^ that the maximum 
downward shift of K^'^^ occurs at wcu = 1- For larger 
wcu, the value of K^'^^{wcu) increases, showing smaU 
oscillations, and saturates for wqu > 6 at a value below, 
though close to, the UIP-MCA constant K^'^^ for the 
undecorated step. A similar tendency is also observed 
experimentally for Cu coverage less than 1.2 ML (cf. Ref. 
I^Ip) but the magnitude of experimental UIP-MA shift, 
jAi^ul, is always larger than ^ \AKy^{wcu = 1)|- 



IV. CONCLUSIONS 

The present tight-binding calculations performed for 
a single step on an ultrathin (001) fee Co slab with in- 
plane magnetization show that the magnetic anisotropy 
of vicinal Co films depends on various factors. In par- 
ticular, the inclusion of crystal field splittings has been 
proved to be vital for the theoretical results to reproduce 
well the experimentally observed decrease of the uniax- 
ial in-plane anisotropy energy induced by decorating the 
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Co steps with Cu atoms. Thus, the generally assumed 
explanation that the shift of the uniaxial anisotropy en- 
ergy is due to the local hybridization between the Co 
and Cu atoms, is confirmed. However, this happens not 
through the extra available hopping between Co and Cu 
atoms, but rather because the potential coming from the 
decorating Cu atoms removes the splitting which exists 
between the on-site energies of d orbitals, with differ- 
ent in-plane orientation, on Co atoms at the edge of the 
undecorated step. The fact that the crystal field split- 
tings affect mainly the anisotropy contribution on the Co 
step-edge atom points to the important role of horizontal, 
in-plane, bonds. Similar conclusion has previously been 
drawn from the experimental results through a general 
analysis based on the Neel's pair-bonding model. 

The investigated shift of the uniaxial anisotropy in- 
cludes also a significant contribution from the local relax- 
ation of atomic structure aromid the step. This finding is 



in a gree ment with the results of earlier studies (cf., e.g., 
Ref. 123) where strain and structure relaxation in the flat 
films have been shown to have great impact on magnetic 
anisotropy. 

These conclusions remain valid also if the charge trans- 
fer between Cu and Co atoms around the step edge is 
present since its effect on the uniaxial anisotropy energy 
of the vicinal Co films has been shown to be very weak. 
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